201801.00632V1 


chinaXiv 


“1476.: 


doi 10.3969/j.issn.1673-4254.2017.11.08 


ChinaXiv 合 作 期 刊 


J South Med Univ 2017, 37(11): 1476-1483 


基础 研究 


血管 内 超声 斑点 的 概率 模型 建立 及 应 用 
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摘要 :超声 交点 是 由 人 体 组 织 中 散射 体 的 反射 信号 相干 作用 所 形成 的 ,其 概率 分 布 与 生物 组 织 的 结构 信息 密切 相关 , 即 不 同 的 
组 织 结构 所 产生 斑点 的 概率 分 布 形式 不 同 。 根 据 血管 内 超声 (IVUS ) 斑 点 的 概率 分 布 特性 ,本 文 提 出 用 斑点 的 伽 马 混合 模型 和 
高 斯 混合 模型 分 别 描述 血管 内 超声 组 织 中 的 钙化 斑 块 软 斑 块 和 正常 血管 区 域 。 通 过 KS 检验 ,KL 散 度 和 相关 系数 等 指标 分 


析 ,发 现 钙化 斑 块 和 正常 血管 区 域 的 斑点 概率 分 布 符合 高 


斯 混 


合 模型 ,而 软 斑 块 更 接近 伽 马 混合 模型 。 在 此 研究 基础 上 ,本 文 


提出 一 种 结合 邻 域 信息 的 概率 混合 模型 ,用 于 IVUS 图 像 斑 块 分 和 


声 影 响 减少 。 
关键 词 :斑点 ;血管 内 超声 ; 斑 块 ;混合 模型 ; 邻 域 信息 


1, 与 现 有 的 概率 混合 模型 比较 ,分 割 精度 大 大 地 提高 , 且 受 噪 


A probability model for analyzing speckles in intravascular ultrasound images to 


facilitate image segmentation 
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Abstract: Ultrasonic image speckles result from the interference of the reflected signals by the scatters in the detected tissue. 


The physical characteristics of the speckles are closely correlated with the structures of the biological tissues, and the 
probability distribution of these speckles differs across different tissues. Based on the probability characteristics of 
intravascular ultrasound (IVUS) speckles, a Gamma mixture model and Gaussian mixture model are proposed to describe the 


calcified plaque, soft plaque and normal vascular regions on IVUS images. Using KS test, KL divergence and correlation 
coefficient analysis, we found that the probability distributions of the speckles generated by calcified plaques and normal 


blood vessels were better described by the Gaussian mixture model while the speckles caused by soft plaques were described 
better by the Gamma mixture model. Based on this finding, we propose a probability mixture model combining neighborhood 


information for plaque segmentation on IVUS images. Compared with the existing probabilistic mixture model, the 
segmentation accuracy was greatly improved with a reduced noise. 
Keywords: speckles; intravascular ultrasound; plaque; probability mixture model; neighborhood information 


超声 成 像 是 利用 超声 波 在 组 织 传播 过 程 中 的 回 波 
信号 进行 成 像 。 在 超声 波 传播 过 程 中 ,由 于 人 体 软 组 织 
声 阻抗 的 不 均匀 性 和 空间 分 布 的 随机 性 ,形成 了 大 量 随 机 


型 部 分 成 熟 型 和 部 分 解决 型 站。 完全 成 熟 型 是 超声 成 
像 中 最 常见 的 斑点 模型 ,反映 组 织 或 者 成 像 区 域 包含 大 
量 散 射 粒子 ,其 结构 小 于 或 等 于 超声 波长 。 


分 布 的 散射 粒子 ,其 散射 波 在 传播 过 程 中 由 于 相干 干涉 
和 相 消 干涉 而 形成 超声 图 像 中 亮 暗 不 定 的 斑点 模式 "。 
斑点 统计 分 布 特性 与 人 体 脏 器 组 织 中 散射 粒子 的 数量 、 
散射 强度 等 密切 相关 ,与 人 体 脏 器 组 织 解剖 结构 特征 直 
接 关联 。 故 斑点 的 物理 特征 反映 组 织 的 结构 特征 。 依 
据 人 体 组 织 中 每 分 辩 单 元 内 散射 粒子 数量 和 粒子 散射 
性 质 ,超声 斑点 分 为 4 种 类 型 :完全 成 熟 型 .完全 解决 
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目前 存在 许多 对 斑点 特征 进行 描述 与 应 用 的 方法 ， 
主要 有 基于 统计 "4 .基于 模型 "信号 处 理 汪 ”和 几何 
方法 上 7 。 其 中 ,基于 统计 的 方法 因 其 简单 有 效 性 逐渐 
引起 人 们 重视 ,并 得 到 了 广泛 的 发 展 。 人 们 已 从 射频 信 
号 (RF) 和 包 络 信号 (超声 图 像 ) 两 个 方面 研究 斑点 统计 
分 布 模型 ,通过 统计 模型 的 参数 估计 ,包括 最 大 似 然 估 
计 、 最 大 后 验 估计 (MAP) 和 EM 算法 中 ,获取 各 器 官 组 
织 特 征 。 斑 点 服从 一 定 概 率 分 布 ,IVUS 斑 块 所 产生 的 
斑点 属于 完全 成 熟 型 斑点 ,基于 完全 成 熟 型 斑点 ,学 者 
们 提出 了 3 种 模型 : 瑞 利 混合 分 布 模型 (RMM) Naka- 
gami 混合 分 布 模型 (NMM) 和 伽 马 混合 分 布 模型 
(GMM)。 其 中 文献 [19] 用 RMM 模 型 描述 斑点 并 对 油 
脂 斑 块 .纤维 斑 块 和 钙化 斑 块 进行 分 类 ,首先 ,用 EM 算 
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法 估计 斑点 的 概率 模型 参数 ,得 到 概率 混合 模型 , 然 > 
法 信 计 玫 的 概率 模 数 得 到 概率 混合 模型 ,然后 p(x18)= Snipi(xl0) ™ 
再 结合 斑 块 组 织 的 纹理 特征 对 斑 块 分 类 。 文 献 [20] 提 和 


出 Nakagami 混 合 分 布 模型 ,该 模型 是 瑞 利 混合 模型 的 
推广 : 瑞 利 分 布 仅 有 一 个 自由 参数 ,而 Nakagami 分 布 有 
2 个 自由 参数 。 这 两 个 模型 根据 超声 射频 信和 号 形成 的 物 
理 过 程 推导 得 出 ,应 用 到 临床 B 型 超声 图 像 分 析 , 其 效 
果 并 不 如 另 一 种 模型 ,这 就 是 文献 [2] 中 提 到 的 伽 马 混 
合 模型 ,该 文中 分 别 用 伽 马 混 合 模型 . 瑞 利 混合 模型 和 
Nakagami 混 合 模型 来 估计 IVUS 图 像 中 组 织 的 概率 分 
布 , 采 用 KL 和 KS 检验 评价 实验 结果 ,发现 伽 马 混合 模 
型 优 于 Nakagami 混 合 模型 ,Nakagami 混 合 模型 优 于 瑞 
利 混合 模型 。 

上 述 研究 认为 斑点 服从 单一 类 型 分 布 ,并 日 指出 用 
伽 马 混合 模型 撒 述 斑点 概率 分 布 最 好 ,本 文 在 其 研究 基 
础 上 ,采用 高 斯 混合 模型 和 伽 马 混合 模型 这 2 种 分 布 类 
型 ,分别 对 IVUS 网 像 中 斑 块 与 正常 血管 组 织 进行 模型 
估计 ,在 不 同 概率 成 分 个 数 取 值 情况 下 ,分 别 用 2 种 概 
率 混合 模型 对 斑 块 和 正常 血管 图 像 进行 概率 拟 合 ; 然 
后 ,用 KS KL 检验 和 相关 系数 等 三 种 指标 评价 实验 结 
果 ,发现 :钙化 斑 块 和 正常 血管 区 域 的 斑点 概率 分 布 符 
合 高 斯 混合 模型 ,而 软 斑 块 更 接近 伽 马 混合 模型 。 说 明 
针对 不 同 组 织 , 可 分 别 选 择 高 斯 混合 模型 和 伽 马 混合 模 
型 对 其 进行 描述 。 本 文 利 用 所 研究 的 概率 模型 ,构建 高 
斯 伽 马 混合 模型 ,并 提出 一 种 结合 邻 域 信 息 的 概率 混合 
模型 ,采用 最 大 后 验 概率 准则 ,实现 ITVUS 图 像 中 各 组 织 
图 像 分 割 ,并 与 文献 [21-23 | 中 混合 模型 用 于 图 像 分 割 
的 方法 比较 ,证 明 本 文 方法 的 优越 性 。 


1 方法 与 流程 

实验 方法 包含 3 个 部 分 :首先 简要 描述 本 文 所 采用 
的 概率 混合 模型 ;然后 介绍 对 混合 模型 效果 进行 评价 的 
3 种 指标 ;最 后 说 明 结合 邻 域 信息 的 混合 模型 用 于 图 像 
分 割 的 具体 方法 。 实 验 流程 如 图 1 所 示 。 
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1 实验 流程 图 


Fig.1 Experimental flow chart. 


1.1 概率 混合 模型 
混合 模型 是 指 随机 变量 * 的 概率 密度 表示 为 : 


其 中 > me=1,9=(r mv0 0v) ,表示 混合 
模型 的 参数 集合 。 即 混合 模型 由 W 个 分 量 组 成 ,每 个 
分 量 的 权重 系数 为 7; 。 知 每 一 个 分 量 都 满足 高 斯 分 
布 , 则 称 其 为 高 斯 混合 模型 ,表达 式 为 : 


p(xl9)= 之 二 os 9 (2) 

参数 集 @= (zo,.….,Ty0,.….,0%) ,Ow =(Ww0n) o 

若 每 一 个 分 量 都 满足 伽 马 分 布 , 则 称 其 为 伽 马 混合 
模型 ,表达 式 为 : 

一 > 

Pr(xl9) 之 BT(a) 

X=0,0,B>0 (3) 

其 中 T(x) 为 欧 拉 伽 马 函 数 , 定 义 : 
T(x)= I exp(-1)qdt , x>0。 参 数 集 

@=(m mxg gw) ,Ou =(ay,B) o 

在 这 里 ,已 知 样本 数据 和 概率 密度 表达 形式 ,利用 
EM 算法 可 以 估计 出 未 知 参数 @ 的 值 ,从 而 得 到 感 兴 
样本 的 混合 分 布 模型 。 
1.2 实验 效果 评价 指标 

本 文采 用 3 种 指标 对 概率 模型 估计 的 效果 进行 评 
佑 ,分 别 是 :KS 检验 KL 散 度 值 和 相关 系数 值 ”。 

KS 检 验 体现 累积 概率 密度 间 的 差异 ,表达 式 为 : 

Due( 忆 ,FJ)=sSup F(x,) F(x) (4) 


其 中 ,五 表示 实际 累积 概率 密度 值 , 屎 . 是 模型 估 
计 的 累积 概率 密度 值 ,sup 是 上 确 界 , 该 式 反 映 实际 累积 
概率 与 模型 估计 累积 概率 间 差 异 的 最 大 值 。KS 检 验 可 
以 判断 实验 样本 是 否 满足 某 种 分 布 ,由 数理 统计 知识 可 
知 : 当 Dj 夺 7TH 时 ,认为 样本 服从 某 种 分 布 是 合理 的 ,本 
文 取 显著 性 水 平 为 0.9, 当 样本 个 数 小 于 100 时 , 查 表 得 到 
阔 值 TF ; 当 样 本 个 数 大 于 100 时 , D, 6,=1.23/Nn ,7 为 
样本 个 数 。 只 要 Due 入 Du。, 即 可 认为 所 提出 的 概率 模 
型 是 合理 的 。 在 KS 检 验 中 : De 越 小 ,概率 模型 估计 效 
果 越 好 。 

KL 散 度 体现 实际 概率 密度 值 与 模型 估计 概率 密度 
值 之 间 的 差异 ,表达 式 为 : 

_~& Pp A(X) 

Dia(pisfi) Zp,(x)log pe) (5) 

其 中 p, 表示 实际 概率 密度 值 ,f 表示 模型 估计 概 
率 密度 值 。 Du 表示 的 是 概率 密度 间 的 差异 ,因此 , Du 
越 小 ,概率 模型 估计 效果 越 好 。 

相关 系数 体现 实际 概率 密度 值 与 模型 估计 概率 密 
度 值 间 相 关 性 的 大 小 ,在 以 上 KS KL 检验 的 基础 上 , 相 


exp(-x/D) ， 
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关系 数值 越 大 ,概率 模型 估计 效果 越 好 。 表 达 式 为 : 


AA 


N N 


DE 


实验 数据 的 实际 概率 密度 由 高 斯 核 平滑 概率 密度 
估计 函数 * 求 得 , 它 是 直方 图 概率 密度 估计 的 拓展 ,已 
被 广泛 应 用 于 数据 概率 密度 估计 ,在 matlab 中 可 以 直接 
调用 依据 该 原理 所 编写 的 ksdensity 函数 对 数据 进行 概 
率 密度 估计 。 

本 文通 过 以 上 3 种 指标 对 概率 混合 模型 估计 的 结 
果 进 行 评判 , 当 某 一 种 概率 模型 估计 结果 的 3 种 指标 中 
2 个 或 以 上 都 是 相应 指标 中 最 优 值 时 ,说明 这 种 模型 的 
估计 效果 更 好 。 

1.3 图 像 分 割 模 型 

传统 混合 模型 由 单一 的 分 布 类 型 混合 得 到 ,本 文 借 
鉴 文献 [25] 中 混合 模型 的 创建 方法 ,用 高 斯 和 伽 马 2 种 
分 布 类 型 构建 高 斯 伽 马 混合 模型 ,模型 表示 如 下 : 

nf (x1O,), 1 <k<m 
Perle@u- 国人 mk<M 

其 中 ,f(x) 表示 x 服从 高 斯 分 布 , (x) 表示 x 服 
从 伽 马 分 布 。 

传统 概率 混合 模型 的 图 像 分 割 方法 没有 考虑 像 
素 间 的 空间 关系 ,所 以 可 能 出 现 以 下 情况 (图 2)。 


(7) 


xD =1 | x(1,D =1 | x1,3)=1 


x(3,D =1 | x(3,2) =1 | x(3,3) =1 


x2.D =1 


x(2,2) =0 


图 2 噪声 图 
Fig.2 Noisy image. 


如 图 2 所 示 ,A.B 图 的 中 心 点 均 为 噪声 点 ,用 传统 
的 概率 混合 模型 会 将 其 与 周围 像素 分 为 不 同 的 类 ,为 消 
除 噪 声 点 对 图 像 分 割 结果 的 影响 ,文献 [22] 中 提 到 对 像 
素 点 后 验 概 率 添加 邻 域 信息 ,但 改善 效果 不 大 ;文献 
[23] 中 提出 噪声 平滑 因子 的 思想 ,但 会 出 现 平滑 过 渡 ， 
细节 丢失 等 问题 。 考 虑 到 后 验 概率 与 像素 先 验 概 率 mr 
和 条 件 概 率 p( xl@,) 直接 相关 ,本 文 结合 像素 间 空 间 关 
系 对 混合 模型 作 以 下 改进 : 

(1) 对 像素 先 验 概率 x, 添加 空间 信息 。 例 如 ,在 图 
2 中 ,A 图 里 所 有 像素 本 应 属于 同一 类 ,但 因 噪声 点 
x(2,2 ) 的 影响 ,可 能 将 其 与 周围 像素 分 为 2 类 。 加 入 
空间 信息 后 ,所 有 像素 点 都 受 周围 像素 影响 ,消除 噪声 
点 干扰 ,所 有 像素 分 为 同一 类 。 


(2) 对 条 件 概率 pfxl@ ) 添加 空间 信息 。 依 据 传统 
概率 混合 模型 ,相同 灰 度 值 的 像素 具有 相同 的 条 件 概 
率 。 在 图 2 中 ,图 A 中 心 像素 与 图 B 周 围 像 素 灰 度 值 相 
同 ,但 它们 不 属于 同一 类 ,传统 混合 模型 不 能 正确 区 分 
两 者 。 加 入 空间 信息 后 , 除 灰 度 值 外 ,所 有 像素 还 受 周 
二 像素 约束 ,从 而 能 更 好 地 对 像素 进行 划分 。 具 体 实 现 
方法 如 下 : 


pxleJ= DR p(x,) (8) 
meN, i 
Dz, 
Nj = (9) 
DD 002 
=1 9， 


k=1lmed, 


其 中 z, 是 EM 算法 中 的 中 间 变 量 ,此 处 z, 计 算 如 下 : 


0 下 四 (10) 
2 DR Pe) 

N 表示 第 i 个 像素 的 邻 域 , 9 = N -全 ,表示 第 i 个 
像素 的 邻 域 但 不 包含 第 i 个 像素 。 w, 是 邻 域 像素 的 权 
重 , R, 是 权重 的 归 一 化 系数 : 

R= Yo, (11) 

本 文采 用 3 x 3 的 均值 模板 对 像素 添加 空间 信息 约 
束 , w, 计算 如 下 : 


1 dnr 
ed (9 


0- 人 (13) 


4 表示 邻 域 像素 六 与 中 心 像素 ;之 间 的 距离 ,由 公 
式 (12) 可 以 看 出 , 邻 域 像素 距 中 心 像素 越 远 ,其 权重 越 
小 , 越 近 则 权重 越 大 。M.N 分 别 表示 模 板 窗口 的 长 和 
宽 , 本 文 窗口 大 小 为 3x3 , 即 MxN=9 , 故 56=2。 可 计 
算出 均值 模板 (图 3)。 


0.155 | 0.176 | 0.155 
0.176 | 0.199 | 0.176 


四 p( 


0.155 | 0.176 | 0.155 


图 3 均值 模板 
Fig.3 Mean template. 


为 避免 无 噪声 图 像 应 用 均值 滤波 时 造成 图 像 部 分 
细节 丢失 问题 ,在 本 文中 ,运用 均值 滤波 时 ,增加 噪声 判 
断 。 当 中 心 像 素 灰 度 值 是 邻 域 像素 中 最 大 或 最 小 值 时 ， 
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认定 它 是 噪声 点 ,使 用 均值 滤波 ;否则 ,不 对 中 心 像素 点 
作 任 何 改变 。 这 样 既 能 很 好 地 结合 像素 点 空间 关系 ,有 
效 抑制 噪声 干扰 ,又 能 保持 图 像 细 节 , 从 而 提高 图 像 分 
割 精 度 。 

图 像 的 分 割 标准 是 最 大 后 验 准 则 。 首 先 ,用 高 斯 伽 
马 混合 模型 对 TVUS 图 像 进行 估计 ,获得 每 一 个 像素 点 
的 后 验 概率 ,然后 根据 最 小 错误 率 贝 叶 斯 决策 完成 每 个 


块 区 域 ,分 别 计算 出 两 种 混合 模型 的 3 种 指标 值 , 最 终 
对 12 幅 图 像 实验 结果 取 平 均值 。 比 较 两 种 概率 模型 的 
实验 结果 ,选择 高 斯 伽 马 混合 模型 ,采用 结合 邻 域 信息 
的 混合 模型 的 方法 ,利用 最 大 后 验 概率 估计 准则 完成 
IVUS 组 织 分 割 。 
2.1 IVUS 各 组 织 概率 混合 模型 分 析 

在 本 节 实 验 过 程 中 ,首先 ,选择 IVUS 图 像 感 兴趣 


像素 的 分 类 ,从 而 完成 图 像 的 初步 分 割 。 因 初步 分 割 图 
像 中 仍然 会 存在 一 些 未 被 正确 分 割 的 孤立 点 ,在 此 基础 
上 ,采用 天 近邻 法 对 分 制图 像 作 后 处 理 ,经 过 后 处 理 ,图 
像 分 割 结果 会 有 较 大 改善 。 


2 结果 

本 文 分 别 用 高 斯 混合 模型 , 伽 马 混合 模型 对 IVU 
图 像 中 的 组 织 进行 概率 模型 估计 。 实 验 数 据 取 自 临 床 
数据 库 ,在 临床 医生 指导 下 ,从 实验 网 像 中 选取 感 兴趣 
区 域 , 感 兴趣 区 域 包括 钙化 斑 块 软 斑 块 和 正常 血管 ,其 
中 钙化 斑 块 选 自 12 幅 IVUS 图 像 , 软 斑 块 选 自 10 幅 
IVUS 图 像 ,正常 血管 选 自 8 幅 IVUS 图 像 。 为 了 分 析 不 
同 混合 模型 和 概率 成 分 个 数 对 实验 结果 的 影响 ,采用 
KL .KS 检验 和 相关 系数 评价 模型 效果 , 即 对 各 类 组 织 
感 兴趣 区 域 构建 混合 概率 模型 ,计算 3 种 指标 值 ,再 分 
别 取 平 均值 。 比 如 :从 12 幅 TVUS 图像 中 ,选取 钙化 斑 


CD 


区 域 ,包括 钙化 斑 块 、 软 斑 块 和 正常 血管 。 根 据 斑 块 回 
声 的 强 弱 ,以 血管 壁 外 膜 回声 为 标准 , 斑 块 定义 如 下 

(1) 软 斑 块 : 软 斑 块 的 回声 弱 于 血管 壁 外 膜 的 回 
声 。 包 括 脂 质 . 疲 松 和 细胞 性 纤维 组 织 。 其 中 脂 质 常常 
表现 为 周围 有 高 回声 区 或 者 纤维 帽 覆 盖 的 低 回声 或 者 无 
回声 区 ,而 踊 松 和 细胞 性 纤维 组 织 则 表现 为 低 回声 区 。 

(2) 钙 化 斑 块 :钙化 斑 块 的 回声 强 于 血管 壁 外 膜 回 
声 ,而 且 后 方 有 清楚 的 声 影 , 钙 化 组 织 占 管 腔 周 径 大 于 
90°, 

图 4 为 IVUS 实 验 图 像 。 通 过 EM 算法 ,分 别 用 两 
种 混合 模型 对 IVUS 各 组 织 估计 公式 (2) (3) 中 未 知 参 
数 ,获得 模型 估计 概率 估计 密度 值 f. 。 在 这 里 ,由 2.2 
中 所 介绍 的 ksdensity 函数 求 出 各 组 织 实际 概率 密度 值 
p, ,结合 公式 (4) (5) (6) ,计算 人 S KL 和 相关 系数 值 ， 
得 到 两 种 概率 模型 在 各 组 织 概 率 估计 的 定量 分 析 结 
(图 5)。 


图 4 实验 原 图 (红色 区 域 为 实验 感 兴趣 区 域 ) 
Fig.4 The original image showing the area of interest (encircled in red). A-C: IVUS of calcified plaques, soft Plaques 
and normal vascular in Cartesian coordinates and polar coordinates, respectively. 


在 图 5 中 ,每 一 行 表示 不 同 组 织 两 种 概率 混合 模型 
的 同一 评价 指标 值 ,在 第 一 行 KS 值 中 ,可 以 看 到 : 当 M> 
2 时 ,各 组 织 成 分 的 KS 值 均 小 于 相应 闷 值 ,符合 概率 模 


型 假设 ;在 M=2 时 ,钙化 斑 块 和 软 斑 块 混合 模型 KS 值 
比 阔 值 大 , 因 IVUS 图 像 中 斑 块 区 域 面 积 较 小 , 即 试验 中 
样本 数 较 少 , 故 概率 拟 合 精度 不 够 。 
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图 5 各 评判 指标 对 比 图 (每 一 列 由 上 到 下 分 别 对 应 各 组 织 相应 指标 图 ) 


Fig.5 Line chart for comparing each evaluation index of different tissues. 


在 图 5 中 ,每 一 列表 示 同 一 组 织 两 种 概率 混合 模型 
的 不 同 评价 指标 值 ,在 A 列 所 对 应 钙化 斑 块 各 指标 值 
中 ,比较 高 斯 混合 模型 和 伽 马 混合 模型 三 种 指标 值 , 可 
以 看 到 ,高 斯 混合 模型 KS KL 值 都 比 伽 马 混合 模型 小 ， 
同时 相关 系数 值 都 比 伽 马 混合 模型 大 ,这 表明 在 IVUS 
钙化 斑 块 组 织 中 ,高 斯 混合 模型 优 于 全 马 混合 模型 , 即 
高 斯 混合 模型 更 符合 钙化 斑 块 。 以 此 类 推 ,从 B.C 列 图 
可 得 出 相应 结论 : 伽 马 混合 模型 更 符合 软 斑 块 , 高 斯 混 
合 模 型 更 符合 正常 血管 。 

在 上 述 对 概率 混合 模型 选择 分 析 基 础 上 ,进一步 说 
明 概率 成 分 个 数 M 值 对 概率 佑 计 的 影响 。 考 察 图 5 中 
每 一 幅 折 线 图 ,可 知 : 随 着 M 值 增 大 ,KS 值 和 KL 值 旦 下 
降 趋势 ,而 相关 系数 值 也 呈 上 升 趋势 。 这 表明 :概率 成 
分 个 数 M 值 越 大 ,概率 拟 合 的 效果 越 好 ;但 当 M 取 值 在 
5 附近 时 ,曲线 变化 趋 于 平缓 ,这 说 明 M 取 值 为 5 时 混合 
模型 的 估计 效果 接近 理想 。 因 此 实际 应 用 中 为 了 减 小 


算法 的 运算 时 间 , 模 型 个 数 可 以 在 5 附近 进行 选取 。 
2.2 结合 邻 域 信息 的 概率 混合 模型 用 于 图 像 分 割 

在 这 里 ,本 文 提出 结合 邻 域 信息 的 混合 模型 用 于 图 
像 分 割 ,以 仿真 图 、 自 然 图 像 科 IVUS 图 像 为 实验 对 象 。 
对 于 仿真 图 和 自然 图 像 ,本 文 用 结合 邻 域 信息 的 高 斯 混 
合 模型 完成 图 像 分 割 ,并 与 现 有 混合 模型 方法 ”比较 ， 
证 明了 本 文 方法 有 效 性 (图 6.7)。 对 于 IVUS 图 像 ,分 
别 用 现 有 混合 模型 ”结合 邻 域 信息 的 高 斯 混合 模型 
和 高 斯 伽 马 混合 模型 完成 图 像 分 割 ,并 给 出 不 同 组 织 成 
分 的 概率 映射 图 ,概率 映射 图 反映 了 分 制图 像 区 域 属于 
某 一 种 组 织 的 概率 ,能 帮助 医师 更 直观 地 看 出 图 像 分 割 | 
的 各 组 织 结构 (图 8.9)。 


3 讨论 
针对 以 上 图 像 分 割 的 结果 ,通过 定性 .定量 分 析 , 比 
较 本 文 分 割 方法 与 现 有 分 割 方法 2 2 的 优 劣 性 。 其 中 ， 
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图 6 仿真 图 分 割 结果 

Fig.6 。 Segmentation results of a simulated image. A-B: 
Simulation image and the experimental image with speckle 
noise; C-F: Segmentation results with Gaussian mixture model 


in 3 previous studies22 and in this study respectively. 


图 7 自然 图 像 分 割 结 果 
Fig.7 Segmentation results of natural image. A: Original image; 
B-E: Segmentation results with Gaussian mixture model in 3 


[21-23] 


previous studies ” and in this study respectively. 


图 8 钙化 斑 块 IVUS 分 割 结 果 
Fig.8 Segmentation results of an IVUS image with calcified plaques. A: Original image; B-F: Segmentation results with 
Gaussian mixture model in previous studies 2 and with Gaussian mixture model and Gaussian-Gamma mixture 
model in this study; G-I: Probabilistic mapping of the vascular lumen, calcified plaques and shadow region. 


定量 分 析 中 利用 到 结构 相似 性 ` 误 分 率 和 峰值 信 品 比 这 
3 种 指标 。 

在 图 6 中 ,图 C 和 图 D 抗 噪 性 较 差 ,图 像 轮 廓 受 只 
声 腐 蚀 严重 。 图 已 和 图 保持 了 图 像 轮 廓 ,但 图 了 中 仍 
残留 较 多 噪声 ,而 图 F 则 有 效 地 抑制 了 噪声 干扰 。 计 算 
现 有 方法 和 本 文 方法 实验 结果 的 结构 相似 性 \ 误 分 率 和 
峰值 信 噪 比 ( 表 1)。 

可 知 ,文献 [23] 中 方法 和 本 文 方法 结果 远 好 于 前 两 
种 方法 ,对 比 文献 [23] 和 本 文 方法 ,本 文 方法 在 结构 相 
似 性 和 峰值 信 品 比 上 比 文献 L23] 中 方法 提高 了 分 别 为 


2.2% 和 29.4%, 此 外 , 误 分 率 比 该 方法 降低 了 68.0%, 可 
见 本 文 方法 比 现 有 方法 有 更 大 优势 。 

在 图 7 中 ,图 B 和 图 C 受 背景 中 噪点 影响 较 大 ,B 图 
中 相机 细节 丢失 ,人 物 轮廓 含 较 多 噪点 。C 图 中 人 物 和 
相机 细节 分 割 效果 比 B 好 ,但 后 方 建筑 物 没有 与 周围 背 
景区 分 开 。D 图 整体 受 周围 噪声 影响 较 小 ,但 细节 丢失 
较 多 ,比如 和 人物 轮廓 有 毛刺 ,照相 机 细节 模糊 ,相机 支架 
与 背景 融合 ,而 E 图 则 弥补 了 D 图 的 缺陷 。 计 算 现 有 方 
法 和 本 文 方法 实验 结果 的 结构 相似 性 如 下 表 2。 

比较 表 中 各 方法 分 割 结果 的 结构 相似 性 ,本 文 方法 
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图 9 软 斑 块 IVUS 分 割 结 果 


Fig.9 Segmentation results of another IVUS image with soft plaques. A: Original image; B-F: Segmentation results with 
[21-23] 


Gaussian mixture model in Previous studies 


and with Gaussian mixture model and Gaussian-Gamma mixture 


model in this study; G-I: Probabilistic mapping of the vascular lumen, soft plaques and shadow region. 


表 1 各 方法 定量 评估 


Tab.1 Quantitative evaluation of each method 


Solution Structural Misclassification Signal- 
Ea similarity index rate noise ratio 
Method in paper2 0.8793 0.3052 14.0122 
Method in paper'™ 0.7673 0.3927 10.6389 
Method in paper'™ 0.9727 0.0603 21.1914 
My method 0.9943 0.0193 27.4144 


表 2 各 方法 定量 评估 


Tab.2 Quantitative evaluation of each method 


Solution Structural similarity index 
Method in paper2 0.7853 
Method in paper.™ 0.7555 
Method in paper'™ 0.7833 
My method 0.8150 


比 文献 [21] 中 方法 提升 了 3.8%, 比 文献 [23] 中 方法 提 
升 了 4.1%。 无 论 视觉 效果 还 是 理论 数据 分 析 , 本 文 方 
法 均 优 于 现 有 方法 。 

在 图 8 中 ,实验 原 图 中 有 3 类 成 分 :血管 腔 \ 钙 化 斑 
块 和 阴影 区 域 。 观 察 到 ,D、E、F 图 分 割 结果 明显 优 于 
BC 图 ,其 中 ,D 图 中 血管 腔 区 域 噪声 点 较 多 ,而 E 下 图 
则 干净 很 多 ,说 明 本 文 结合 邻 域 信息 的 混合 模型 法 的 分 


制 结果 明显 优 于 现 有 混合 模型 法 ;E 图 和 F 图 都 能 将 钙 
化 广 块 很 好 地 分 割 出 来 ,但 F 图 中 蝶 声 点 比 EE 图 少 很 多 ， 
区 域 分 割 更 完整 ,说 明 高 斯 伽 马 混合 模型 优 于 高 斯 混合 
模型 。G HI 图 分 别 为 血管 腔 、 钙 化 谢 块 和 阴影 区 域 所 
对 应 的 概率 映射 图 。 以 图 G 举 例 说 明 ,图 中 ,红色 区 域 
表示 血管 腔 , 其 概率 值 最 大 , 蓝 色 部 分 为 背景 区 域 ,概率 
值 较 小 。 

在 图 9 中 ,实验 原 图 中 有 3 类 成 分 :血管 腔 、 软 斑 块 
和 阴影 区 域 。 可 以 看 到 ,B 图 没有 将 阴影 区 域 从 血管 腔 
区 域 分 离 干 净 , 且 噪点 较 多 ,C 图 误 将 血管 内 腔 区 域 划 
分 为 阴影 区 域 ,D.E 图 较 好 地 分 割 开 各 区 域 ,但 两 图 中 
血管 内 腔 区 域 噪声 点 较 多 ,区 域 分 割 不 完整 ;而 F 图 则 
很 好 地 弥补 了 这 些 不 足 , 分 割 区 域 连贯 且 干 净 ,进一步 
说 明 本 文 高 斯 伽 马 混合 模型 法 优 于 现 有 混合 模型 法 。 
GHJI 图 分 别 为 血管 腔 、 软 斑 块 和 阴影 区 域 所 对 应 的 概 
率 映 射 图 。 以 图 五 举例 说 明 ,图 中 ,绿色 线 所 标记 的 红 
色 区 域 为 软 斑 块 ,其 概率 值 最 大 , 浅 蓝 色 区 域 为 血管 腔 ， 
深蓝 色 区 域 为 阴影 区 域 ,两 者 概率 值 都 较 小 。 

血管 内 超声 斑点 的 像素 强度 符合 一 定 的 概率 分 布 ， 
本 文 用 伽 马 混合 模型 和 高 斯 混合 模型 来 估计 钙化 斑 块 
和 软 斑 块 所 产生 斑点 和 正常 血管 腔 的 概率 分 布 ,KS 检 
验 结果 表明 这 两 种 分 布 能 很 好 地 描述 斑点 和 正常 血管 
腔 的 统计 特性 。 此 外 ,用 KS、KL 和 相关 系数 这 3 种 指标 
评价 两 种 概率 模型 佑 计 结 果 , 得 出 结论 :钙化 斑 块 和 正 
常 血管 腔 用 高 斯 混合 模型 拟 合 其 分 布 效果 更 好 ;而 软 斑 
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块 则 用 伽 马 混合 模型 拟 合 其 分 布 效果 更 好 。 基 于 这 个 
结论 ,在 具体 IVUS 图 像 分 割 中 ,依据 图 像 中 各 组 织 分 布 
情况 ,确定 用 哪 一 种 概率 混合 模型 进行 图 像 分 割 。 在 
IVUS 图 像 分 割 中 ,本 文 用 上 述 两 种 模型 构建 高 斯 伽 马 
混合 模型 ,并 结合 像素 空间 关系 ,以 最 大 后 验 概率 准则 
对 图 像 进行 分 割 。 经 试验 结果 验证 ,本 文 所 提出 的 结合 
邻 域 信息 的 混合 模型 优 于 现 有 混合 模型 ,不 仅 提 高 了 分 
割 精 度 ,而 且 有 效 地 抑制 了 噪声 干扰 。 但 因为 基于 概率 
混合 模型 的 分 割 方法 主要 依据 的 是 图 像 的 灰 度 信息 ,所 
以 易 将 IVUS 图 像 中 纹理 特征 区 别 较 小 , 灰 度 值 接近 的 
区 域 划分 为 同一 类 。 考 虑 到 灰 度 值 接近 但 属于 不 同 组 
织 的 区 域 间 有 分 界 区 域 , 故 解决 方法 有 :(1) ,临床 医生 
可 对 比 原 图 和 分 割 图 像 确定 灰 度 值 相近 区 域 边界 ,从 而 
将 两 区 域 区 分 开 ; (2) 将 本 文 方法 与 其 它 基于 图 像 轮 廓 
的 分 割 方法 结合 , 减 小 对 灰 度 值 的 依赖 ,提高 图 像 分 割 


质量 。 
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